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Higgs boson production rates at liadron colliders are reviewed with particular emphasis on progress 
in the calculation of higher order QCD effects. Emphasis is placed on the uncertainties in the 
predictions for Higgs boson production. A firm understanding of these uncertainties is crucial for 
extracting new physics signals. 



g . I. INTRODUCTION 

o : 

■ The search for the Higgs boson is one of the fundamental goals of the Tevatron and the LHC. At the Tevatron, 
^ [ discovery relies on the observation of the Higgs signal in many channels, since the rate is quite small. In order to 
D . disentangle a Higgs boson from the background, it is crucial to have reliable predictions for both the signal and 

' the background and to well understand the uncertainties in the theoretical predictions. The case is somewhat 
■ . different at the LHC, since the Higgs boson production rate is significantly larger than at the Tevatron. Here 
I ' discovery is more straightforward, and the question becomes what can we learn about the underlying physics. 
This again requires a firm understanding of the production rates and decay patterns. 
' In this paper, we review the predictions for the production of the Higgs boson in the Standard Model at 
J> I the Tevatron and at the LHC. The current status of QCD radiative corrections to Higgs boson production is 
I examined critically in order to make estimates of the uncertainties involved in the predictions. The next-to- 
leading order (NLO) results are now available for all of the dominant production channels and the forefront of 
activity has moved to the calculation of next-to-next-to-leading order (NNLO) corrections and the resummation 
I ' of leading and next-to-leading order logarithmic effects. The implementation of these higher order effects in 
, Monte Carlo programs and the comparison with fixed order perturbative calculations remains an active area of 

■ investigation. 

I ■ II. LHC 

Oh! 

^ I Many studies have been made of the capability of the LHC to observe the Higgs boson in a variety of channels. 
• • . With 30 fb~^, a Standard Model Higgs boson will be observable at the 5a level over the entire mass range, 
100 GeV < Mh < 1 TeV. With 100 the Higgs boson will be observable in at least 2 channels in the 

, same mass range. By combining various channels, some measurements of Higgs couplings will then be possible. 

■ The conclusions of these studies, however, typically rely on the use of the lowest order cross sections only.Q 
Ci I The NLO QCD results exist for all the relevant Higgs production and decay channels, but not for many of the 

backgrounds. 



A. Gluon Fusion 



The dominant production mechanism for the Higgs boson at the LHC is gluon fusion. The NLO QCD 
corrections are well known, both in the Mt — > oo limit and with the inclusion of the complete Alt dependence 
of the result. [H A convenient parameterization of the results is given by the K factor. 

We have explicitely included the dependence on fi to emphasize the fact that the K factor is typically quite 
sensitive to the renormalization/factorization scale fi. The Mt — > oo limit provides an extremely accurate 
description of the full rate at NLO, as can be seen in Fig. By including the exact Mh and Mt dependence in 
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TABLE I: pp -> at ^ = 14 TeV (CTEQ5, set4; fi = Mh) 
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FIG. 1: Complete NLO result for inclusive Higgs boson production at the LHC, \/s — 14 TeV, with a renormaliza- 
tion/factorization scale fi = Mh (solid). The dashed line is obtained by combining the complete lowest order result with 
the K factor computed in the limit Mt — > oo. 



the lowest order result and multiplying by the K factor computed in the Mt — > oo limit, ((Tmt^oo), the resulting 
approximation to the next-to-leading order rate is extremely accurate all the way up to Mh ~ 1 TeV. Table 1 
shows the dependence of the NLO cross section on various input parameters. The exact NLO calculation has a 
small dependence on rUb through the b quark loop. 

The NLO corrections to the LO rate are quite large, increasing the cross section by about a factor of 2. In 
addition, the scale dependence remains significant. The LHS of Fig. 2 compares the complete LO order result 
with the NLO result. The bands represent a variation of the renormalization scale from ^ < fj, < 2Mh- Note 
that there is no overlap between the bands labelled LO and NLO and so, in this case, the variation of the scale 
/Lt appears to be a poor indicator of the uncertainty of the result. 

The accuracy of the NLO K factor computed in the Mt — > oo limit has encouraged two groups to undertake 
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FIG. 2: The curved labelled LO (NLO) on the LHS are the lowest order (next-to-leading-order) result for pp — > /i at the 
LHC with the renormalization scale varied from -^^ < fi < 2Mh- The curve labelled NLO-SVC on the LHS includes 
the terms of Eq. 1 with i — j — 1, along with the NLO virtual corrections. The curve labelled NNLO-SVC on the RHS 
includes the terms of Eq. 1 with i = 1, 2, 3 and of Eq. 2 with j — 3, plus the NNLO virtual corrections. From Ref. [4]. 



the calculation of the NNLO contribution to inclusive Higgs production in this limit. ^ Using an effective 
theory corresponding to an infinite top quark mass, the NNLO virtual corrections reduce to 2-loop Feynman 
diagrams, instead of the 3-loop diagrams they would be in the complete Mt dependent calculation. These virtual 
contributions to the NNLO rate have been computed by Harlander.Q 

At present, the existing NNLO results for the inclusive Higgs production rate are incomplete and make an 
assumption which the authors term the "soft approximation" . This approximation includes the leading terms 
as z = M^/s 1, (where is the gluon-gluon center of mass energy). These leading z ^ 1 contributions are 
of the form, 

^(1-.), -M,3 . (1) 

These terms are expected to provide the bulk of the NNLO corrections. The vaUdity of the soft approximation 
can be tested at NLO by comparison with the complete calculation. Inclusion of the leading soft terms of Eq. 
(1) with i — 1, (the curve labelled NLO-SV of Fig. 2), shows that the soft plus virtual contributions alone 
underestimate the exact NLO result by ~ 15 — 20%. 

In order to obtain a more accurate approximation to the complete rate, the sub-leading coUinear contributions 
can also be included. These terms are of the form 

log^(l-z), J = 1,2,3. (2) 

The leading coUinear contributions at each order (j=l for NLO and j — 3 for NNLO) have been found in 
Refs. [|[ Q. In addition, the sub-leading coUinear contributions (j = 1, 2 at NNLO) can be estimated from the 
resummation calculation of Ref. From Fig. ^, we see that including the coUinear log(l — z) contribution, 
along with the virtual contribution and the soft terms of Eq. 1 (with i — 1), provides an excellent approximation 
to the fuU NLO resuU (the curve labelled NLO-SVC in Fig. 2). 

Using the soft plus coUinear approximation to the NNLO result, i ~ 1, 2, 3 in Eq. (1) and j = 3 in Eq. (2), 
yields the results shown on the RHS of Fig. | (labelled NNLO-SVC). We see that the NNLO corrections are 
large, leading to a i^T factor between 2.5 and 3. The bands correspond to varying the renormalization scale 
between Mh/2 < fj, < 2Mh. The scale dependence is only slightly reduced from that of the NLO result. 

Harlander and Kilgorej^ included also the sub-leading coUinear terms of the form log^(l — z) and log(l — z), 
to obtain the solid curves shown in Fig. ^. The differences between the 3 upper curves in Fig 3 is due taking 
different approximations for the unknown sub-leading coUinear contributions to the NNLO result and can be 
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FIG. 3: (a) NNLO result for the K factor for pp /i at ^ = 2 TeV with = Mh.. (b) NNLO resuh for the K factor 
for pp ^ ft at ^ = 14 TeV with ii = A'h. From Ref [3]. 
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FIG. 4: Uncertainty on the production cross section times branching ratio to NLO at the LHC for various Higgs channels. 
^J, is varied from Mh/2 < ^ < 2Mh,, Qs(Mz) from .115 < as{Mz) < .123 and nib from 4.5 GeV < mt{mt) < 5 GeV. 



interpreted as an estimate of the uncertainty of the result. As is clear from this figure, the coUinear contribution 
from the log^ (1 — z) terms is numerically quite large, (since the dotted curve labelled "soft" omits the log-' (1 — z) 
contributions). The inset in the RHS of Fig. 2 also show the importance of the coUinear contributions. 

Another important issue is the question of NNLO parton distribution functions (pdfs). At present only partial 
NNLO pdfs exist. Catani et.al. use the NNLO pdfs of Ref. 0, while Harlander and Kilgore utilize CTEQ5 
NLO pdfs. The inclusion of NNLO pdfs (instead of NLO pdfs) decreases the rate by roughly 8%.Q Clearly 
a complete NNLO calculation with complete NNLO pdfs is needed before we can begin to extract precision 
results. At present, the best estimate is that there is still approximately a 35% uncertainty in the prediction 
due to scale dependence, unknown NNLO terms, incomplete knowledge of the NNLO pdfs, and our knowledge 
of as- 

An important outcome of higher order calculations is the Higgs boson pT spectrum. At lowest order, the Higgs 
boson is produced with no transverse momentum. At higher orders in ag , the effects of soft and soft plus coUinear 
gluon emission from the initial state partons are numerically significant. At low px, the usual factorization 
approximation fails and large logarithms of the form a" log™(M|/p|n) appear. These large logarithms can be 
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FIG. 5: Comparison of the Higgs pr spectrum derived using soft gluon resummation at low pr, matched with the exact 
calculation to 0(Qs) at high pr (solid), with PYTHIA (dashed). The 3 solid curves are estimates of the uncertainty 
from unknown NNLO contributions. From Ref. [8] 



resummed to give a result which is valid at low pT- At an intermediate value of px, the resummed form can 
be matched with the exact matrix element calculation to ©(a^), valid at large pt- This result is shown in Fig. 
^. The three solid curves in this plot represent an attempt to estimate uncalculated NNLO contributions to 
the resummed result. Ref. Q estimates a ±10% uncertainty due to these unknown terms. When the complete 
NNLO calculation is available, it will be possible to remove much of this uncertainty. 

Since experimental searches rely strongly on Monte Carlo programs, it is important to understand how soft 
gluon emission is included in these programs. Monte Carlo programs typically produce the Higgs px spectrum 
using parton showering, which correctly reproduces the spectrum at low py, but underestimates the rate at 
higher pT, as can be seen in Fig. ^. The correct spectrum at high pT (as determined from the exact matrix 
element calculation) can be obtained using PYTHIA by judiciously adjusting the arbitrary renormalization 
scale. § 

The NLO rate for Higgs plus 1- jet production at the LHC has been computed in the Mt oo limit. |^ 
The corrections increase the rate by a factor of 1.5 — 1.6 and are almost constant over a large range of Mh, 
rapidity, and pT- As with inclusive Higgs production, the renormalization/factorization scale uncertainty is 
significantly reduced by the inclusion of the NLO contributions, although the residual uncertainty is still rather 
large, ^ ±20%. The NLO QCD results for all of the 2- and 3- body Higgs decays have existed for some time and 
are conveniently implemented in the FORTRAN code HDECAY.|^0| Fig. ^ shows the variation of the inclusive 
Higgs production cross section calculated to NLO multiplied by the NLO branching ratios to various channels 
as the renomalization scale, /i, is varied from Mft/2 < /i < 2M;j, as{Mz) is varied within the LEP 1-a limit, 
.115 < as[Mz) < .123 and mf, is varied within the particle data group range, 4.5 GeV < mb{mb) < 5 GeV. The 
dominant source of uncertainty in the results of Fig. 4 is the renormalization/factorization scale dependence. 
By measuring combinations of final states, the uncertainty on the predictions can, however, be significantly 
reduced. |l^ The theoretical uncertainty on the branching ratios alone is considerably smaller than on the 
product crB.flll 
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FIG. 6: (a) pp —> 2 jets + h from the weak boson fusion sub-process at ^/s = 14 TeV (dashed) and background from 
the gg — > ggh sub-process calculated exactly (solid) and in the Mt oo limit (dotted), (b) Same as (a), but with cuts 
designed to enhance the weak boson fusion contribution. From Ref [14]. 

B. Vector Boson Fusion 

Vector boson fusion can be used to measure the WWh and ZZh couplings at the LHC.[0 The NLO QCD 
corrections to pp -^2- jets +h through the vector boson sub-process are quite smaU, and the uncertainly 
on the production rate is estimated from the small scale dependence to be ~ 1 — 2%. 

It is necessary to separate the signal, qq qqh ( which probes the WWh and ZZh couplings), from the 
background, gg — > ggh (which depends only on the tih Yukawa coupling) and first enters at l-loop.Q Since 
the dominant contribution to the background arises from gluons in the initial state, the background is enhanced 
by the large gluon luminosity at the LHC. The gg ggh contribution to the Higgs plus 2- jet signal has been 
computed both in the Mt — > oo limit and retaining the full Mt mass dependence and the results are shown in 
Fig. ^. The rrit ^ oo limit is a good approximation to the full result only for M^ < 2Mt. The weak boson 
fusion process produces well separated, forward jets, while the jets from the gluon fusion sub-process are more 
isotropic. Using cuts designed to enhance the vector boson fusion contribution, shown in Fig 6(b), it is clear 
that it will be possible to separate the gluon fusion sub-process from the vector boson fusion contribution for 
Mh < 2Mt. The weak boson fusion processes dominate over the gluon initiated processes by about a factor of 
3 to 1 after applying the appropriate cuts and for Mh < 2Mt. 

C. Associated Production 

A Higgs boson in the mass range around 120 — 140 GeV is particularly difficult to observe at the LHC since 
the preferred channel, h — > 77, suffers from a small rate and large backgrounds. In this region, the associated 
production channel, pp tth, may be useful to confirm an elusive Higgs signal. Although the production rates 
are small, ~ .5 — .8 pb, the signature with the final state W^W~bbbb is spectacular. This process is of particular 
interest since it can be used to measure the tth Yukawa coupling. 

The tth process proceeds predominantly through gluon fusion at the LHC. The complete NLO results have 
been found in Ref. jl^ and are shown in Fig. [3. The NLO predictions show a significantly reduced scale 
dependence and increase the rate by roughly 20% from the LO predictions over the entire intermediate Higgs 
mass range. 

III. TEVATRON 

Higgs boson production rates at the Tevatron are much smaller than at the LHC, but with the increased 
luminosity of Run II, it may be possible to observe a Higgs signal for a Higgs mass below around 180 GeFjie). 
The dominant production mode is gluon fusion, with a cross section between 1.0 and 0.2 pb at y/s — 2 TeV for M^ 
in the 120—180 GeV region. Gluon fusion, however, suffers from large QCD backgrounds to the dominant h ^ bb 
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FIG. 7: (a) Dependence of crLO,NLo(pp —> tth) on the renormalization/factorization scale fi, at y/s = 14 TeV, for 
Mh = 120 GeV. (b) onlo and Uio for pp ith as a function of Mh, at y's = 2 TeV, for ^ = Mt and = 2Mi. From 
Ref [15]. 



decay channel. The gluon fusion production mechanism may, however, be useful for 140 < Mh < 180 GeV, 
when combined with the h WW* decay channel. 

The most likely discovery channel at the Tevatron for Mh < 140 GeV is the associated production of Wh or 
Zh, where an efficient trigger is provided by the leptonic decay modes of the vector bosons. For Higgs bosons 
in the Mh ^ 120 GeV region, the associated production with a top quark pair may also be observable. The 
status of the NLO QCD corrections to these processes is briefly discussed below and the rates at the Tevatron, 
including NLO QCD corrections for all channels, are shown in Fig. |[ (The rate for Zh production is about a 
factor of 2 below that for Wh production). 



A. Gluon Fusion 



At the Tevatron, gluon fusion contributes roughly 65% of the total Higgs production cross section for 
120 GeV < Mh < 180 GeV. The NLO corrections to inclusive Higgs production are large and positive for 
all values of the Higgs mass. The NNLO corrections to gg ^ h at the Tevatron have been computed by Har- 
lander and KilgoreQ and by Catani ei.a/.Q in the soft plus coUinear approximation described above and are 
shown in Fig. ^. Given the large numerical value of these partial corrections, a complete NNLO calculation is 
essential before reliable predictions can be made in this channel. 



B. Associated Production, pp Wh, Zh 



The Wh and Zh channels are the most promising discovery channels at the Tevatron for Mh < 140 GeV. 
The NLO rate for pp Wh is shown in Fig. || and is around .1 — .2 pb. (This figure does not include the W and 
h decay branching ratios.) The NLO QCD corrections are the same as those for Drell-Yan and increase the rate 
by about 30% from the lowest order prediction. The dependence of the NLO corrected rate on the choice of 
parton distribution functions is quite small, but there remains about a 12% uncertainty in the prediction due 
to the residual renormalization/factorization scale dependence. 

Since the dominant decay of a Higgs boson below Mh ~ 140 GeV is to bb pairs, the irreducible background 
processes to pp — > Wh and pp — > Zh are pp — > Wbb and pp — > Zbb. These background processes have been 
calculated to NLO in Rcf. p9|| and the results implemented in the Monte Carlo program, MCFM. The NLO 
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FIG. 8: NLO rates for Higgs production in pp collisions at ^/s — 2 TeV, evaluated at the renormalization/factorization 
scale fi = Mh- (The process pp tih is evaluated at /i = 2mt). 
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FIG. 9: The m^s distribution of the W^{—* VeC^) bh background to the Higgs signal, {pp Wh), at LO and NLO at 
^/s = 2 TeV . The fe-quark tagging efficiency is not included. From Ref. [19]. 



corrections are large and positive and change the shape of the bb distribution near the peak, as can be seen in 
Fig. ^. The K factors for the background Wbb and Zbb processes are larger than those for the Wh and Zh 
signals and have not been included in the studies of Ref. ||l9[ . 



C. Associated Production, pp tth 

At the Tevatron, the associated production of tth proceeds primarily through qq annihilation. While the rate 
for tth production is small, the signature is distinctive. Unlike at the LHC, at the Tevatron the invariant mass 
distributions of the final state bb pairs from the tih signal have rather a different shapes from the background 



9 




FIG. 10: (a) Dependence of aLO,NLo{pp — * tth) on the renormalization/factorization scale fj,, at -^s = 2 TeV, for 
Mh = 120 GeV. (b) onlo and crio for pp tth as a function of Mh, at = 2 TeV^ for /i = Mt and ^ = 2Mt. From 
Ref [20]. 



and preliminary studies suggest that it may be possible to observe this channel at the Tevatron.|lj 

The next to leading order results have been computed recently by two groups, with excellent agreement. 
The NLO result shows a reduced scale dependence from the lowest order result and a slightly reduced cross 
section from the lowest order prediction. For example, for Mh — 120 GeV and /x = Mt, the NLO total cross 
section is reduced to 4.86 ± 0.03 fh from the lowest order prediction of 6.868 ± .002 fh. The reduction is much 
less dramatic at = 2Aft, as can be seen from Fig. Only for renormalization/factorization scales larger 
around than /i = 2Mf + Mh is the NLO cross section larger than the lowest order rate. Combining the residual 
scale dependence with the error from the parton distribution functions (~ 6%) and from (7%), we estimate 
the uncertainty on the theoretical prediction as about 12%. 



IV. CONCLUSIONS 



NLO corrections to all Higgs production channels of interest at hadron colliders have now been completed. 
These NLO predictions show a significantly reduced renormalization/factorization scale dependence from the 
LO predictions, leading to increased confidence in the validity of the predictions. Complete NNLO corrected 
rates for inclusive Higgs production should be available soon. Consistent NNLO calculations will, however, 
require structure functions derived to NNLO, which are not yet available. 
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